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S. Verga, A. Knigavko, and F. Marsiglio 
Department of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2J1 

Recent energy dispersion measurements in several families of the hole-doped copper oxides have 
revealed a kink in the energy vs. momentum relation. These have tentatively been identified as 
due to electron phonon coupling. We invert this data directly to determine the bosonic spectral 
function; the kink gives rise to a singular function in the phonon energy region. 



The determination of mechanism for superconductivity 
in the high temperature oxides has occupied researchers 
for the past fifteen years. The most definitive signature 
for determining mechanism in conventional superconduc- 
tors traditionally has been the measurement of the sin- 
gle particle tunneling I-V characteristic^, and the con- 
comitant inversion procedure^. Researchers^ have re- 
ported some success with this procedure for the high 
temperature superconductors; nonetheless the applica- 
bility of such a procedure is unknown outside a weak 
coupling electron phonon framework, and the process 
is complicated in the superconducting state due to the 
non-isotropic nature of the order parameter. Other pro- 
cedures have been suggested, such as the inversion of 
the normal state optical conductivity*, the conductivity 
in the superconducting state in conjunction with neu- 
tron scattering data^, and the inversion of photoemission 
data^ in the superconducting state. Recent very high 
resolution photoemission measurements on a variety of 
cuprate materials have suggested that a sizeable electron 
phonon coupling exists^, and the possibility of inverting 
this data (in the normal state) has been re-opened. 
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FIG. 1. Dispersion of occupied states vs. momentum for 
LSCO with three different doping concentrations. The insert 
shows a schematic which puts this data in the context of a 
parabolic band. A modification due to electron phonon cou- 
pling occurs near the Fermi surface. The region corresponding 
to the data of the main figure is highlighted in bold on the 
left side of the parabola. 



In this paper we outline an inversion procedure, deter- 
mine the 'bosonic spectrum' to which the electrons are 
coupled, and assess the possibility that these bosons are 
the phonons. The data for the dispersion for three dif- 
ferent dopant concentrations in LSCO is reproduced in 
Fig. 1^. Particularly for the underdoped sample there 
is a well-defined kink which occurs at approximately 70 
meV. Lanzara et alJ attributed this kink to an electron 
self energy effect due to coupling to phonons. We wish to 
investigate this claim based on some microscopic models. 

We first examine the result obtained from 'standard' 
phonon models, namely the Einstein and Debye mod- 
els. Each in turn is used to model the electron phonon 
spectral function, a^Fiy), and then, within the standard 
framework^'^, the electron self energy = -|- 

iYj2{'jj) is obtained (at T = 0): 



(1) 



from which the electron dispersion Ek can be obtained: 



Ek = ek + ^i{Ek), 



(2) 



where is the bare (with respect to electron phonon 
interactions) quasiparticle energy. The model spectrum, 
along with the real part of the self energy and the dis- 
persion is plotted in Fig. 2(3) for the Einstein (Debye) 
models, respectively. Note that for some values of €k the 
dispersion as shown is multivalued. This is because mul- 
tiple poles exist in these regions. In fact, the spectral 
function generally evolves in the following manner* (we 
use the Einstein model depicted in Fig. 2 for simplicity): 
at low energies a pole exists just above the real axis (i.e. 
with infinitesmal width). What is not depicted in Fig. 
2 or 3 is that the weight of this pole (i.e. the residue) 
goes to zero as \ek\ — > oo- So, beyond about |efc| « 140 
meV, the weight of this pole becomes very small. The 
actual energy {Ek) at this point becomes very nearly the 
Einstein frequency. In a sense the electrons and phonons 
have 'hybridized' and this branch, which started out very 
'electron-like' near ii'fe « is now very 'phonon-like'. The 
branch just below -80 meV in Fig. 2c is more or less ir- 
relevant, since the self energy has a very large imaginary 
part (see Fig. 2b). Finally, as \ek\ continues to increase 
the lowest branch begins to dominate. In this limit this 
branch becomes very 'electron-like' albeit with a finite 
width. 
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FIG. 2. The (a) boson spectrum vs. frequency, (b) real 
and imaginary parts of the electron self energy vs. energy 
below the Fermi level, and (c) the resulting dispersion vs. 
bare quasiparticle energy. These figures are for the Einstein 
model for phonons, with uje = 80 meV and \e = 0.5. See 
text for a discussion of the multivalued portion. 



S(a; + i6) = dv a^F{v) -2iTi (n{u) + 
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(3) 



where il^{x) is the digamma function and n{i') is the Bose 
distribution function. Similarly, for the Dcbyc model 
there is a well pronounced shoulder, but still no clear 
kink. To determine what sort of spectrum does lead to a 
kink we utilize the procedure described below. 
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FIG. 4. Temperature dependence of (a) the real part of the 
self energy, and (b) the dispersion relation for the Einstein 
model shown in Fig. 2. Note that the singularity is washed 
out by finite temperature. 



Debye phonon spectrum (T=OK) 



Self energy 



U- 

^ 0.4 
0.0 



{3.) . 


(Oq = 80 meV 




1-0 = 0.8 








100 200 -100 

(0 (meV) CO (meV) 

Quasiparticle dispersion 



200 



> 

E-100 



-200 





Ek 


(c) 




Ek 




coq = 80 meV 






Xq = 0.8 






E|,=ek+Re[E{E^)] 











200 



100 

lEkl (meV) 

FIG. 3. Same as in Fig. 2, except for the Debye model, 
with u!d = 80 meV and Xd = 0.8. 



In addition, in the case of the Einstein spectrum there 
is an unphysical singularity. While this is easily washed 
out by temperature, thermal effects still do not reconcile 
theory with experiment. Example calculations are shown 
in Fig. 4. We use the full finite temperature expression 
for the self energy®: 



As a preliminary analysis we note that the following 

function fits the data in the underdoped regime (where 
the kink is particularly prominent) reasonably well: 




if \Ek\ < u)D 
if \Ek\ > ujd 



(4) 



This parametrization is motivated by the following con- 
siderations: we wanted an analytical form that would 
contain an explicit 'kink', we wanted as simple a form 
as possible, and we wanted a form with which we could 
analytically perform the Kramers-Kronig integration to 
obtain the imaginary part of the self energy. We have 
had to use an additional fit for the high energy region, 
to relate to k — kp- For simplicity, we have used a 
linear fit, i.e. ek = xk' , with k' = \k — fci?|. Inclusion of a 
quadratic term (see Ref.^" for a discussion of quadratic 
corrections) results in a correction which is negligible. 
The proportionality constant x has units of eV A and is 
related to the bare Fermi velocity [x = hvF)- For the un- 
derdoped sample, the proportionality constant has been 
determined by the requirement that, at high energy, the 
full dispersion becomes the bare one (the model self en- 
ergy decreases to zero). Because the kink is washed out 
by increased doping, using the same function to fit the 
data in the optimally doped and especially in the over- 
doped regime is difficult. To overcome this problem we 
renormalize the momentum such that the three experi- 
mental curves overlap in the high energy range. We then 
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use the fit for the undcrdopcd data to calculate the pro- 
portionality constant x for the high energy fit for the 
optimally doped and overdoped sample. 

The result of using three fitting parameters, A and u)n, 
and the background slope x is shown in Fig. 5; it is 
clear that the fit is very good. The important feature, 
captured by the fit, is the initial linear dependence of 
energy on wavevector, followed by an abrupt change at 
some characteristic frequency. In Fig. 5b we show the 
extracted real and imaginary parts of the self energy from 
the data. We show the model fit, the actual data once the 
high energy part is extracted, and the smoothened curve 
used to carry out the Kramers-Kronig analysis. We find 
parameter values of A = 0.89, ujd = 72 meV, and vp = 
4.1 eV-A/h (= 6.2X10'' cm/s). Varying these parameters 
'by hand' results in a less than 10 % change in A, for 
example, before the fit becomes visually poor; we thus 
regard this as a rough measure of the uncertainty. 

Quasiparticle dispersion versus momentum 
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FIG. 5. (a)The fit (solid curve) to the measured disper- 
sion data (shown with symbols) for the underdoped case. We 
have used Eq. (3) with ujd = 72 meV and A = 0.89. The 
dotted line shows the linear fit that results using x = 4.1 eV 
A. (b) Real part of the self energy vs. energy, as extracted 
from the data (symbols), the model fit (solid curve), and the 
smoothened fit to the data (dashed curve). The insert shows 
the imaginary part obtained through Kramers-Kronig analy- 
sis. 



The imaginary part of the self-energy can be deter- 
mined from the real part through a Kramers-Kronig in- 
tegral: 

E2{lo + iS) - £2(00 + i5) = - 



d.'^^. (5) 



Substituting Eq. (4) results in 



2 / uj 

^2(0; + i6) - ^2(00 + iS) = -A LUnf — 



(6) 



where f(x) 
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With the standard 



approximations^'^^, one can relate the imaginary part of 
the self energy to the underlying electron phonon spec- 
trum. The required result is: 



IT uM 



(7) 



The result from the fit shown in Fig. 5 is plotted in 
Fig. 6. Note that f{x) is the same function which ap- 
pears in the Hartree-Fock calculation for the free electron 
gas^^, and as is well known, its derivative has a logarth- 
mic singularity at the characteristic frequency, as shown. 
Aside from the logarithmic singularity, the spectrum is 
peculiar (as an electron phonon spectrum) in that it is 
linear at low frequencies, and has a long tail at high fre- 
quency. A model-independent calculation (i.e. without 
the fit given in Eq. 4) for each of the doping concentra- 
tions provided in Fig. 1 is also shown in Fig. 6. There 
is clearly a rounding of the singularity^^, although the 
overall strength of the interaction, as indicated by the 
area under the spectrum, is of the same magnitude. 
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FIG. 6. The bosonic spectral functions which result from 
an inversion of the data for underdoped, optimally doped, 
and overdoped samples (solid, short-dashed, and long-dashed 
curves, respectively). Parameters for the fits are given in 
Table. 1. The model based on Eq. 3, designed specifically 
for the underdoped case, 5 — 0.07, is also shown (dot-dashed 
curve), along with that expected for spin fluctuations^'*. In 
the spin fluctuation model of^* there is considerable spectral 
weight at high frequencies, including weight up to 400 meV 
(not shown). Note that in the inverted experimental results 
there have slight negative portions, which are set to zero in 
the subsequent analysis in the text. 
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Several questions arise from the results of Fig. 6. First, 
what sort of coupling strength do these spectra represent, 
as measured by a superconducting critical temperature, 
assuming that these are utilized in an Eliashberg-type 
analysis ? To answer this question, one would ideally 
like to perform a calculation with an order parameter 
with d-wave symmetry. However, for this to be mean- 
ingful one would require a series of results as shown in 
Fig. 1 for various directions in the Brillouin zone. Then 
the same inversion scheme would result in a momentum 
dependent spectrum, which would then be used in a d- 
wave Eliashbcrg equation. Since this information is lack- 
ing, we simply use an s-wavc calculation, with the direct 
Coulomb interaction set to zero. This would apply if 
the repulsion was primarily short-range; then the d-wave 
symmetry would be unaffected by it. On the other hand, 
there is no guarantee that the momentum dependence 
of the data acquired in this way would necessarily lead 
to d-wave superconductivity. We proceed in this sim- 
plistic way nonetheless, and obtain the results for the 
three doping concentrations summarized in Table 1. As 
judged both by the value of and the parameter A, 
the coupling strength decreases as the doping increases. 
Probably data from more dopant concentrations is re- 
quired before one can assess how significant this trend 
is, and what additional physics may causing To to de- 
crease in the underdoped regime. We should emphasize 
that the main part of the analysis is done in the normal 
state. Since momentum dependent data is lacking, we 
have tacitly assumed that the coupling strength is inde- 
pendent of momentum; an important refinement will be 
to eventually try to extract this momentum dependence 
and see how it correlates with d-wave pairing. 



dopant level 


5 = 0.07 


5 = 0.15 


5 = 0.22 


A 


0.94 


0.82 


0.51 


a;in(meV) 


42 


40 


38 


Tc (K) 


58 


46 


19 



TABLE I. Spectral function parameters as a function of 

dopant concentration. The parameters A and cuin arc deter- 
mined by numerical integrals^ for the spectral functions shown 
in Fig. 6. All indicators show a somewhat enhanced coupling 
at optimal doping. 



Another interesting question is whether these spectral 
functions represent phonons or not. As pointed out in 
Rcf.^, the most compelling evidence favouring phonons 
is that the frequency domain is consistent with that ob- 
served in neutron scattering. Fig. 6 does show high 
frequency spectral weight, however, and one can ask 
whether these high frequency tails (clearly beyond the 
phonon energies in these materials) rule out phonons as 
a possibility. To address this question we cut off the 
spectrum at 100 meV, and then readjusted other param- 
eters in the fit (by hand) to recover an improved fit to 
the data originally presented in Fig. 1. While the fit 
is never as good as the original one, we find that it is 
sufficiently good to be a plausible possibility. Thus, un- 
fortunately, we are unable to say anything very definite 
on this issue. It is true that spin fluctuations (one of 
the competing alternatives to phonons) are expected to 
have significant spectral weight at higher frequency. An 
example is shown in Fig. 6, taken from Ref.""^^, which has 
considerable spectral weight extending up to 400 meV. 
It is clear that significant spectral weight exists at fre- 
quencies much higher than indicated by the experimen- 
tal data; furthermore, this particular spectrum appears 
to be considerably softer below the 50 meV region than 
the data indicates. Once again more definite statements 
could be possible once a detailed momentum dependence 
of the spectral functions is available. 

Finally, we have focussed on the real part of the dis- 
persion. One might well ask why we didn't examine the 
imaginary part directly. The partial answer to this ques- 
tion is that there arc other (i.e. non-pairing) scattering 
processes which affect the imaginary part of the self en- 
ergy and not the real part (e.g. impurities — see below). 
Even more critically, as was attempted in the original 
angle-resolved photoemission spectroscopy (ARPES) in- 
version work^, one might try to invert the entire spectral 
function. The difficulty is exemplified in Fig. 7, where 
we show an energy distribution curve taken a,t kp for 
the underdoped sample. We have also plotted the en- 
ergy distribution calculated with our model fit for the 
self energy, Eqs. (4) and (6), additionally including an 
energy resolution fimction^"'' and elastic scattering from 
impurities (which affects the imaginary part of the self 
energy, but not the real part). The fit is excellent at low 
energies, but there is a clear and very large discrepancy 
at high energies. The origin of this discrepancy is simply 
not understood at present. It may or may not represent 
new physics (there may be a lack of understanding in the 
analysis of ARPES), but it clearly will not be understood 
with the approach adopted in this paper, which focuses 
on the quasiparticle peak as an essential feature. Mainly 
for this reason we felt that an examination of the real 
part of the self energy (i.e. the dispersion) was the best 
procedure for extracting a potential pairing interaction. 
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FIG. 7. The energy distribution function (spectral function 
times the fermi function convoluted with an energy resolution 
function^^ vs. energy, for the underdoped sample, at kp- The 
impurity scattering rate used was 1/t = 80 meV, and the en- 
ergy broadening was 8 meV. The fit is excellent at low energy, 
but clearly fails at high energy. 

In summary we have used the electron dispersion, as 
measured by ARPES, to extract an electronic self energy 
(real part). Through Kramcrs-Kronig we are able to ex- 
tract the imaginary part of the self energy, from which 
an inelastic scattering spectral function is extracted in 
a straightforward manner. The result is summarized in 
Fig. 6 for a variety of doping concentrations in LSCO. 
The result is clearly compatible with phonons; the ex- 
tracted coupling strength would then be able to account 
for the superconductivity. However, based on the avail- 
able data we are unable to rule out, for example, spin 
fluctuations as a possibility. Measurements over more 
doping concentrations and in different directions in the 
Brillouin zone would aid considerably in narrowing down 
the possibilities. 
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